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In  this  work,  an  energy-harvesting  system  which  extracts  heat  from  an  automotive  exhaust  pipe  and 
turns  the  heat  into  electricity  by  using  thermoelectric  power  generators  (TEGs)  was  investigated.  The 
influences  of  the  number  and  the  coverage  rate  on  the  heat-exchanger  of  the  TEGs  were  explored  via 
simulations.  It  was  found  that  implementing  more  TE  couples  does  not  necessarily  generate  more  power 
in  total,  and  most  of  all  the  average  power  per  TE  couple  decreases  rapidly.  It  is  because  the  wall  tem¬ 
perature  of  the  exhaust  pipe  drops  quickly  along  the  streamwise  direction  and  also  because  the 
downstream  TEGs  contend  for  heat  with  the  upstream  TEGs,  causing  a  reduction  in  the  temperature 
difference  between  the  hot  and  cold  sides  of  the  upstream  TEGs.  Furthermore,  it  was  also  found  that  for  a 
given  total  number  of  TE  couples,  it  is  better  to  retain  a  portion  of  the  heat  exchanger  uncovered  with  TE 
couples  at  the  downstream  side  so  that  the  downstream  wall  of  the  exhaust  pipe  uncovered  with  TE 
couples  becomes  even  hotter  than  the  upstream  wall  covered  with  TE  couples.  Heat  is  consequently 
conducted  from  the  downstream  wall  to  the  upstream  wall  and  successively  to  the  attached  TEGs;  a 
larger  total  power  can  be  thus  obtained. 

©  2013  Elsevier  Masson  SAS.  All  rights  reserved. 


1.  Introduction 

Because  of  the  global  energy  crisis  and  the  environmental  pro¬ 
tection  issues,  energy  recovery  techniques  have  become  signifi¬ 
cantly  demanding  for  a  long  time.  Some  examples  of  energy 
recovery  techniques  are  water  heat  recycling,  heat  recovery  venti¬ 
lation,  heat  recovery  steam  generators,  and  so  on  [1].  Waste  heat 
recovery  by  using  thermoelectric  power  generators  (TEGs)  is 
another  attempt.  TEGs  can  directly  convert  thermal  energy  to 
electrical  energy  and  have  the  advantages  of  light  weight,  no  noise, 
and  no  mechanical  vibration.  Owing  to  these  merits,  TEGs  have 
found  its  potential  in  many  applications,  such  as  space  applications, 
thermal  energy  sensors,  textiles,  etc.  [2—5]. 

Waste  heat  originates  from  various  heat  sources  at  different 
temperatures,  such  as  human  bodies,  plants,  automobiles,  and  so 
on  [5].  Human  body  provides  low- temperature  heat  for  TEGs.  The 
electrical  energy  converted  from  the  human-body  heat  can  be 
applied  directly  to  wearable  electronic  products  such  as  watches 
and  electronic  medical  instruments  [6,7].  These  small-scale  appli¬ 
cations  can  be  operated  by  small  temperature  differences  [8]. 
Waste  heats  generated  from  plants  and  automobiles  are  at 
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medium-to-high  temperatures  and  can  sustain  much  larger  tem¬ 
perature  differences  across  the  TEGs,  resulting  in  large  power 
generation.  A  500  kW  TEG  system  for  waste  heat  recovery  in  a 
pressurized  water  reactor  power  plant  was  developed  by  Central 
Research  Institute  of  Electric  Power  Industry  (CRIEPI)  [9].  In  addi¬ 
tion  to  the  power  generation,  the  reductions  in  the  fuel  consump¬ 
tion  and  in  the  CO2  emission  resulted  by  installing  TEGs  are  also 
noticeable.  Chen  et  al.  [10]  in  use  of  a  Danish  district  heating  system 
model,  estimated  the  reductions  in  a  system  which  integrates  TEGs 
into  a  combined  heat  and  power  production  (CHP)  plant.  The  re¬ 
sults  provide  useful  information  for  developing  energy  policies  and 
estimating  their  environmental  impacts. 

Waste  heat  from  automotive  vehicles  is  considerable  as  well.  For 
a  typical  gasoline-fueled-internal-combustion-engine  vehicle, 
about  40%  of  the  fuel  energy  is  discharged  from  the  exhaust  pipe 
and  about  30%  is  lost  into  the  coolant.  Making  good  use  of  these 
waste  heats  improves  the  energy  efficiency  and  saves  money.  Yang 
[11]  indicated  that  the  consumer  fuel  savings  over  a  three-year 
period  is  about  $400  for  a  23.5  mpg  vehicle,  under  the  assump¬ 
tion  of  $2/gallon,  15,000  miles/yr,  and  a  desired  10%  fuel-economy 
improvement  (the  overall  objective  raised  by  the  US  Department  of 
Energy  in  2004).  The  commonly  utilized  components  in  a  vehicle 
for  implementing  the  TEGs  are  the  radiators  and  the  exhaust  sys¬ 
tem.  Hsiao  et  al.  [12]  established  mathematic  models  and  per¬ 
formed  experiments;  they  found  a  better  performance  can  be 
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Nomenclature 

Qd 

the  total  amount  of  heat  transferred  from  the 

downstream  side 

Cp 

the  heat  capacity  of  air 

Qh 

the  heat  extracted  from  the  hot  side  of  a  TEG 

he 

the  characteristic  convection  heat  transfer  coefficient 

Q.in 

the  input  extractable  energy 

of  the  heat  sink 

Rm 

the  measured  electric  resistance  of  a  TE  couple 

h 

the  convection  heat  transfer  coefficient  for  flat-plate 

Rpn 

the  electric  resistance  of  a  TE  couple 

flow 

ReLc 

the  Reynolds  number  based  on  the  characteristic 

I 

the  electric  current 

length  Lc 

I< 

the  thermal  conductance  of  the  TEG  cuboid 

AT 

the  temperature  difference  between  the  hot  and  cold 

K 

the  thermal  conductivity  of  air 

sides  of  a  TEG 

km 

the  effective  thermal  conductivity  of  a  TEG  chip 

Ta 

the  ambient  temperature 

L 

the  length  of  the  heat  exchanger 

Tc.avg 

the  averaged  temperature  at  the  cold  side  of  a  TEG 

Ljeg 

the  length  of  a  TEG  cuboid 

Tj-i.avg 

the  averaged  temperature  at  the  hot  side  of  a  TEG 

Lc 

the  total  length  of  the  simulated  system 

Tin 

the  inlet  temperature 

M 

the  weight  of  the  TEG  system 

ATieg 

the  temperature  difference  across  a  TE  leg 

m 

the  mass  flow  rate 

V 

the  speed  of  the  vehicle 

N 

the  total  number  of  TE  couples  in  the  system 

w 

the  side  length  of  the  hexagonal  pipe 

Nx 

the  number  of  TE  couples  along  the  x-direction  in  one 

^saved 

the  saved  engine  shaft  power 

TEG  cuboid 

the  measured  Seebeck  coefficient  of  a  TE  couple 

Nz 

the  number  of  the  TE  couples  along  the  z  direction  in 

Wpn 

the  Seebeck  coefficient  of  a  TE  couple 

one  TEG  cuboid 

?7 alternator 

the  alternator  efficiency 

Nu 

the  average  Nusselt  number 

Vo 

the  driveline  transmission  efficiency 

P 

the  power  generation  rate 

^exchanger 

the  efficiency  of  the  exchanger 

^max 

the  maximum  power  generation  rate 

VpCU 

the  Power  Conditioning  Unit  (PCU)  efficiency 

Pr 

the  Prandtl  number 

V  recovery 

the  recovery  efficiency 

Q.  absorb 

the  total  amount  of  heat  transferred  to  the  TEGs 

Vteg 

the  efficiency  of  the  TEG 

Qc 

the  heat  dissipated  at  the  cold  side  of  a  TEG 

Mr 

the  rolling  resistance  coefficient 

obtained  by  attaching  TEGs  to  the  exhaust  system  than  to  the  ra¬ 
diators.  Directly  attaching  TEGs  to  the  original  equipments  in  the 
exhaust  system,  such  as  the  catalytic  converter  [13],  the  muffler,  the 
exhaust  pipe,  and  so  forth,  is  one  way  to  go.  Alternatively  one  can 
construct  a  heat  exchanger  to  extract  heat  from  the  exhaust  pipe  to 
the  integrated  TEGs. 

Depending  on  the  kind  of  the  cooling  system,  the  TE  materials, 
the  system  size,  and  the  operating  conditions,  the  power  genera¬ 
tion  rate  differs  in  the  literature.  The  Hi-Z  Technology,  Inc.  [14] 
attached  TEGs  on  the  outer  surface  of  an  octagonal  cast-steel  heat 
exchanger  with  a  hollow  displacement  body  in  the  center  and  used 
eight  aluminum  water-cooled-heat-sink  assemblies.  By  installing 
swirl  fins  on  the  displacement  body  and  using  seventy-two  HZ-13 
Bi2Te3  modules,  they  obtained  1  kW  from  a  Cummins  NTC  350 
Diesel  engine  operated  at  300  HP  and  1700  RPM.  In  2001,  in  order  to 
solve  mechanical  problems  during  the  road  test,  HZ-14  modules 
were  built.  A  maximum  electric  power  of  900  W  was  obtained  from 
a  Cummins  335  Diesel  engine  operated  at  290  HP  and  2100  RPM 
[15].  Crane  and  Lagrandeur  [16]  built  a  single-layer  high-tempera- 
ture  segmented  TEG  with  the  most  suitable  thermoelectric  mate¬ 
rials  for  various  operating  temperature  ranges.  A  power  as  large  as 
125  W  was  obtained  by  extracting  heat  from  a  hot  air  at  600  °C  and 
at  a  flow  rate  of  45  cfm  and  using  a  cooling  water  at  25  °C  and  at  a 
flow  rate  of  10  1pm.  Karri  et  al.  [17]  simulated  a  TEG-integrated 
sports  utility  vehicle  which  includes  sixteen  HZQW  (quantum- 
well  materials,  Si/SiGe  and  B4/B9C)  TEG  chips.  The  water  coolant 
system  in  the  vehicle  was  also  used  to  cool  the  TEG  system.  About 
450  W  was  generated  at  a  vehicle  speed  of  112.7  km/h  and  about 
1.5%  fuel  savings  was  obtained  after  the  consideration  of  the 
coolant-pumping  power  loss,  blow-down  power  loss,  and  the  po¬ 
wer  loss  to  transport  the  weight  of  the  TEG  system;  the  fuel  saving 
of  the  car  using  Bi2Te3  TEGs  is  negative  nonetheless.  Hsu  et  al.  [18] 
proposed  a  heat  exchanger  mounted  with  eight  Bi2Te3-TEG  chips 
and  employed  eight  air-cooled-heat-sink  assembles.  A  maximum 
power  of  44  W  was  obtained.  In  2012,  they  enlarged  the  heat 


exchanger  with  a  similar  interior  structure  and  added  a  slopping 
block  in  the  inlet  [19].  A  maximum  power  of  12  W  was  obtained 
instead  when  twenty-four  TEG  chips  and  twenty-four  air-cooled- 
heat-sink  assembles  were  implemented.  Their  experimental  results 
surprisingly  suggest  that  a  small-size  heat  exchanger  with  less  TEGs 
can  generate  more  power. 

Most  of  the  studies  mentioned  above  focused  on  the  design  of 
the  heat  exchangers  and  the  development  of  the  TEG  modules. 
However,  according  to  Hsu’s  study,  the  size  of  the  heat  exchanger 
and  the  number  of  the  TE  couples  implemented  are  two  other 
factors  that  may  strongly  affect  the  power  generation  efficiency. 
Therefore,  the  influences  of  these  two  factors  are  of  interest  in  the 
present  work.  The  investigation  was  executed  by  performing  full 
three  dimensional  simulations  in  use  of  the  commercial  software 
ANSYS-Fluent.  A  heat  exchanger  similar  to  that  employed  by  the  Hi- 
Z  Technology,  Inc.  [14]  was  employed.  As  its  total  length  and  the 
coverage  rate  with  the  TE  couples  were  varied,  the  power  genera¬ 
tion  rate  was  computed.  The  optimization  of  these  two  parameters 
is  aimed  at. 

The  rest  of  this  paper  introduces  the  simulation  model  in  Sec¬ 
tion  2,  presents  the  simulation  results  and  discussions  in  Section  3, 
and  finally  gives  the  conclusion  in  Section  4. 

2.  Numerical  experiment 

2.1  Simulation  model 

The  energy  harvesting  system  consists  of  a  heat  exchanger  and 
TEGs  as  shown  in  Fig.  1.  The  heat  exchanger,  composed  of  a  hex¬ 
agonal  pipe,  radial  fins,  and  a  hollow  center  body,  is  connected  to 
the  exhaust  pipe  of  diameter  62  mm  on  both  sides.  The  hexagonal 
pipe  has  an  inscribed  circle  of  diameter  140  mm  and  a  length  of  L.  It 
is  connected  to  the  exhaust  pipe  via  a  divergent/convergent  part. 
The  hollow  center  body,  including  a  circular  pipe  of  diameter 
100  mm  and  length  L  as  well  as  two  bullet-shape  heads  of  length 
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Fig.  1.  The  energy  harvesting  system:  (a)  a  three  dimensional  chart,  (b)  an  x—y  cross- 
sectional  view,  (c)  the  simulation  domain  (one-sixth  of  the  system),  (d)  the  dimensions 
on  the  cross  section  of  x  =  0. 


boundary  condition,  referring  to  Hsu  et  al.  [18],  is  a  uniform  flow  of 
velocity  12  m/s  and  temperature  Tm  =  773  K.  Imposed  at  the  outlet 
boundary  are  zero  gradients  for  velocity  and  temperature.  The  in¬ 
ner  surface  of  the  center  body  is  adiabatic.  The  heat  sink  attached  to 
the  cold  side  of  the  TEG  was  not  simulated  but  modeled  by  a 
characteristic  convection  heat  transfer  coefficient  he.  In  this  study, 
attempted  are  two  values,  600  W/m2  I<  and  1800  W/m2  K,  which 
characterize  the  S2090-40W  heat  sink  produced  by  Alpha  Company 
Ltd.  operated  at  an  air  velocity  of  1  m/s  and  5  m/s  respectively  [21  ]. 
The  rest  of  the  outer  surface  of  the  heat  exchanger  on  the  other 
hand  has  a  convection  heat  transfer  coefficient  of  h  =  60  W/m2  K. 
This  value  is  obtained  based  on  the  following  empirical  equation 
[22]: 


Nu  =  =  Pr1/3 (o. 037 Re°Lf  -  87l)  (1) 

for  flat-plate  flows  with  Lc  =  460  mm  (the  total  length  of  the  energy 
harvesting  system  when  L  =  180  mm;  see  Fig.  Id),  an  air  velocity  of 
30  m/s  (about  the  car  speed  on  the  freeway),  and  air  properties 
evaluated  at  the  ambient  temperature  Ta  =  300K,  where  appear  the 
air  thermal  conductivity  fca,  the  Prandtl  number  Pr,  and  the  Rey¬ 
nolds  number  ReLc.  As  the  heat  exchanger  length  L  varies  from 
50  mm  to  240  mm,  Eq.  (1)  gives  a  value  of  h  ranging  from  53  W/ 
m2  K  to  65  W/m2  K.  The  variation  is  small.  For  convenience,  we 
adopt  a  fixed  value  of  h  =  60  W/m2  K  in  all  the  simulations. 
Finally,  the  exhaust  gas  was  modeled  as  air  whose  properties  vary 
with  temperature,  and  the  TEG  was  modeled  as  a  cuboid  of  thick¬ 
ness  2.8  mm  with  an  effective  thermal  conductivity  of  km  =  3.26  W / 
m  K  [18,19]. 


90  mm,  is  supported  by  the  radial  fins  of  thickness  1  mm,  which  are 
connected  to  the  inner  surface  of  the  hexagonal  pipe  on  the  other 
end.  All  of  the  solid  components  are  made  of  aluminum.  Finally, 
TEGs  of  length  LTeg  (colored  by  red  in  Fig.  1)  are  attached  to  the 
outer  surfaces  of  the  hexagonal  pipe,  starting  from  the  most  up¬ 
stream  edge. 

2.2.  Simulation  method 

Because  of  the  symmetry  of  the  energy-harvesting  system,  only 
one-sixth  of  the  system  needs  simulating.  An  unstructured  mesh 
system,  which  is  composed  of  about  3.5  x  106  cells  and  about 
1.8  x  106  nodes,  was  carefully  built.  The  SIMPLE  algorithm  and 
second-order  upwind  scheme  were  used  to  solve  the  steady-state 
heat  transfer  problem.  The  realizable  k—e  turbulence  model  was 
adopted,  in  which  the  dissipation  rate  equation  is  derived  from  the 
dynamic  equation  of  the  mean-square  vorticity  fluctuation  at  large 
turbulent  Reynolds  number  [20].  The  convergence  criteria  adopted 
herein  are  a  scaled  residual  under  10-3  for  the  momentum  balance, 
10~6  for  the  energy  balance,  and  a  relative  error  under  0.1%  for  the 
total  energy  conservation  of  the  system. 

The  boundary  conditions  are  illustrated  in  Fig.  2.  The  no-slip 
boundary  conditions  are  imposed  at  all  the  solid  walls.  The  inlet 
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2.3.  Power  calculation 

Once  the  simulated  system  has  reached  steady,  the  temperature 
difference  between  the  hot  and  cold  sides  of  the  TEG  cuboid  (TH.avg 
and  Tc.avg)  was  extracted  to  compute  the  power  generation  rate. 
The  generated  power  of  a  TEG  is  generally  expressed  as 

P  =  NapnIATleg  -  l2NRpn,  (2) 

where  ATieg  is  the  temperature  difference  across  the  TE  leg,  N  is  the 
number  of  TE  couples  employed,  I  is  the  electric  current,  apn  is  the 
Seebeck  coefficient,  and  Rpn  is  the  internal  electric  resistance  of  one 
TE  couple.  However,  instead  of  the  temperature  difference  across 
the  TE  leg,  it  is  usually  the  temperature  difference  between  the  two 
ceramic  plates  of  a  TEG  chip  (AT)  that  is  measured.  Besides,  the 
copper  plate,  the  contact  thermal  as  well  as  electric  resistances,  and 
the  heat  loss  due  to  convection  and  radiation  heat  transfers  all  may 
cause  a  loss  in  the  output  power.  We  thus,  in  evaluating  the  power 
generation  rate,  adopt  the  measured  Seebeck  coefficient  am  and  the 
measured  electric  resistance  Rm  (see  Table  1)  of  the  TEG  chips 
employed  by  Hsu  et  al.  [23],  each  of  which  is  composed  of  199 
Bi2Te3-TE  couples. 

On  the  other  hand,  because  the  cross  section  of  a  single  TE 
couple  is  about  3.77  mm  x  3.77  mm,  the  total  number  of  TE  couples 
implemented  are  N  =  6 Nz  x  Nx,  where  Nz  =  LTEG/3.77  mm  and 
Nx  =  W  13.77  mm  with  W  =  48  V3  mm  being  the  side  length  of  the 


Table  1 

The  thermoelectric  properties  of  the  employed 
TEG  chips  [23]. 


a  m 

Rm 


Uniform  air  flow  (z  direction) 
12m/s,  773K 


Fig.  2.  The  boundary  conditions  employed  in  this  study. 


295  |iV/K 
6.48  mQ 
3.26  W/m  K 
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Table  2 

The  simulation  results  when  L 
he  =  1800  W/m2  K. 

=  LjEG  =  180  and 

KATavg 

4836  W 

hiavg 

375.9  K 

Fc.avg 

329.5  K 

ATavg 

46.16  K 

hexagonal  pipe.  The  maximum  power  generation  rate,  which  can 
be  obtained  by  selecting  an  external  load  resistance  equal  to  the 
total  internal  electric  resistance  (. NRm ),  and  the  corresponding 
electric  current  therefore  become 

Pmax  =  3.36  x  10-6-NATa2vg(w)  (3) 

and 

I  =  amATavg/(2  Rm)  (4) 

with  Aravg  being  the  averaged  temperature  difference  between  the 
hot  and  cold  sides  of  the  TEG  cuboid. 

There  is  one  more  issue  that  needs  addressing  before  we  end 
this  section.  In  the  present  simulations,  we  did  not  consider  the 
Joule’s  heat  and  Peltier  heat;  that  is,  an  open-circuit  system  was 
simulated.  The  so-obtained  temperature  difference  ATaVg  was  then 
employed  to  evaluate  the  power  generation  rate  by  using  Eq.  (3).  It 
is  known  however  that  ATavg  actually  varies  with  the  external  load 
resistance  or  the  generated  electric  current  due  to  the  Peltier  and 
Joule  heat.  To  estimate  the  induced  error,  we  take  the  system 
having  L  =  LTeg  =  180  mm  and  he  =  1800  W/m2  I<  for  example.  The 
temperatures  at  the  hot  and  cold  sides  of  the  TEG  and  the  heat 
conducted  to  the  TEG  obtained  from  the  open-circuit  simulation 
are  listed  in  Table  2.  We  now  consider  an  imaged  close-circuit 
system  within  which  the  temperatures  at  the  hot  and  cold  sides 
of  the  TEG  are  as  the  same  as  those  in  the  open-circuit  system.  The 
heat  extracted  from  the  hot  side  (QH)  and  the  heat  dissipated  at  the 
cold  side  (Qc)  of  a  TEG  in  the  close-circuit  system,  namely, 

Qh  =  N«m/TH, avg  +  KA^avg  -  l/2NRm,  (5a) 

Q  c  =  Nam/TCjavg  +  K  ATcWg  +  \ I2NRm  (5b) 

are  found  to  be  Qh  =  5537  W  and  Qc  =  5492  W;  the  values  of  the 
Peltier  and  Joule  heat  at  both  sides  are  thus  701  W  and  656  W 
respectively,  which  are  about  14%  of  the  conduction  heat  /<ATavg 
(I<  is  the  thermal  conductance  of  the  TEG  cuboid). 

3.  Result  and  discussion 

In  the  following  two  subsections,  we  investigate  the  influence  of 
the  heat-exchanger  length  (L)  and  its  coverage  rate  with  the  TE 


Fig.  4.  The  distribution  of  the  temperature  difference  between  the  hot  and  cold  sides 
of  the  TEG  cuboid  in  a  system  in  which  L  =  LTEG  =  180  mm  and  he  =  1800  W/m2  K. 

couples  on  the  generated  power  separately.  In  the  former,  the 
performances  of  the  heat  exchangers  of  various  lengths,  completely 
covered  by  TE  couples,  were  compared.  In  the  latter,  we  fixed  the 
total  number  of  the  TE  couples  and  changed  the  coverage  rate  by 
varying  the  length  of  the  heat  exchanger.  A  comparison  between 
the  completely  and  partially  covered  situations  was  made  in 
Section  3.3.  The  waste-heat  recovery  efficiency  was  discussed  in 
Section  3.4. 

3.1.  Effect  of  the  length  of  the  heat  exchanger  completely  covered 
with  TE  couples 

In  this  subsection,  we  make  the  outer  surfaces  of  the  heat 
exchanger  completely  covered  by  the  TE  couples;  that  is  L  =  LT eg- 
Intuitively  one  expects  the  longer  the  heat  exchanger  length,  the 
more  TE  couples  are  installed,  and  consequently  the  more  power  is 
generated.  Five  different  lengths,  50  mm,  80  mm,  120  mm,  180  mm, 
and  240  mm,  were  attempted.  Fig.  3  shows  the  x—y  and  y—z  cross- 
sectional  temperature  distributions  of  the  system  having 
L  =  lTEG  =  180  mm  and  he  =  1800  W/m2  K.  Fig.  4  shows  the  dis¬ 
tribution  of  the  temperature  difference  between  the  hot  and  cold 
sides  of  the  TEG  cuboid.  As  seen,  because  of  the  cooling  effect  (by 
he),  the  temperature  of  the  exhaust  air  decreases  rapidly  along  the 
streamwise  direction;  so  is  the  temperature  difference  between  the 
hot  and  cold  sides  of  the  TEG  cuboid.  The  temperature  difference 
variation  in  the  x  direction  is  much  smaller  (Fig.  4),  so  we  averaged 
the  temperature  difference  along  the  x  direction  and  observed  its 
variation  in  the  z  direction  as  shown  in  Fig.  5.  It  is  observed  that  the 
temperature  difference  decreases  with  the  downstream  distance 
due  to  the  drop  in  the  fluid  temperature.  Furthermore,  the  tem¬ 
perature  difference  at  a  fixed  downstream  distance  decreases  with 
increasing  heat-exchanger  length.  This  is  because  the  strong  heat 
sinks  mounted  at  the  downstream  TEGs  attract  heat  to  be  con¬ 
ducted  from  the  upstream  side  to  the  downstream  side  via  the  solid 
components  in  the  system.  Consequently  the  longer  the  heat- 
exchanger  length,  the  less  heat  is  conducted  to  the  upstream  TE 
couples  and  the  smaller  the  upstream  temperature  difference  be¬ 
comes.  Finally,  without  surprise  the  stronger  the  heat  sink,  the 
larger  the  temperature  difference  is. 


Lug  =  1 80mm 


Fig.  3.  The  x—y  and  y—z  cross-sectional  temperature  distributions  in  a  system  in  which  l  =  LTEG  =  180  mm  and  he  =  1800  W/m2  K. 
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Fig.  5.  The  distributions  of  the  temperature  difference  between  the  hot  and  cold  sides  of  the  TEG  cuboid  along  the  z-direction  when  the  heat  exchanger  is  completely  covered  with 
TE  couples. 


We  now  calculate  the  total  power  generation  rate  and  the  average 
power  generation  rate  per  TE  couple  according  to  Eq.  (3).  The  results 
are  shown  in  Fig.  6.  As  expected,  the  total  power  increases  with  the 
increasing  heat-exchanger  length  or  the  increasing  number  of  TE 
couples  but  gets  saturated  quickly.  However,  the  average  power  per 
TE  couple  has  an  opposite  variation  trend.  This  is  because  the  total 
power  generation  is  proportional  to  the  total  number  of  TE  couples 
but  to  the  square  of  the  temperature  difference;  increasing  the  heat- 
exchanger  length  increases  the  number  of  TE  couples  but  decreases 
the  temperature  difference  as  discussed  above. 

3.2.  Effect  of  the  length  of  the  heat  exchanger  partially  covered  with 
TE  couples 

We  now  try  to  cover  the  outer  surface  of  the  heat  exchanger  with 
TE  couples  partially.  We  fixed  the  length  of  the  TEG  cuboid,  LTeg»  at 
50  mm  and  varied  the  length  of  the  heat  exchanger.  The  x-averaged 
temperature  differences  between  the  hot  and  cold  sides  of  the  TEG 
cuboid  are  shown  in  Fig.  7.  It  shows  that  the  temperature  difference 
increases  with  the  heat  exchanger  length  and  it  decreases  first  but 


increases  later  near  the  downstream  edge  of  the  TEG  cuboid  along 
the  downstream  distance.  A  further  examination  on  the  wall  tem¬ 
perature  of  the  heat  exchanger  as  shown  in  Fig.  8  finds  the  down¬ 
stream  wall  without  TE  couples  mounted  (black  curves)  has  a  higher 
temperature  than  the  upstream  wall  with  TE  couples  mounted  (red 
curves)  due  to  a  smaller  convection  heat  transfer  coefficient 
(h  =  60  W/m2  I<)  over  there.  Heat  is  therefore  transferred  to  the 
endmost  TEGs  not  only  from  the  hot  fluid  but  also  from  the  hotter 
downstream  wall;  a  larger  temperature  difference  is  consequently 
resulted.  We  distinguished  and  listed  the  amounts  of  heat  trans¬ 
ferred  to  the  TEGs  in  Table  3,  where  Q. absorb  denotes  the  total  amount 
of  heat  transferred  to  the  TEGs  and  Qd  is  the  amount  transferred 
from  the  downstream  side.  As  seen,  Qd  can  reach  as  much  as  13%  of 
Qabsorb  When  he  =  600  W/m2  K  and  17%  when  he  =  1800  W/m2  K. 
Besides,  Q.d  increases  with  the  heat-exchanger  length  but  quickly 
gets  saturated.  This  explains  the  rapid  saturation  observed  in  the 
temperature  difference  and  the  total  power  generation  rate  (Fig.  9). 
From  the  economic  consideration,  an  incomplete  coverage  is  supe¬ 
rior  to  a  complete  coverage  and  a  heat  exchanger  slightly  longer 
than  twice  the  length  of  the  TEG  cuboid  is  sufficient;  the  latter  is 
suggested  from  the  dashed  line  indicated  in  Fig.  9. 


Fig.  6.  The  total  power  (solid  lines)  and  the  average  power  per  TE  couple  (dash  lines) 
against  the  heat-exchanger  length  L  when  the  heat  exchanger  is  completely  covered 
with  TE  couples. 


3.3.  Comparison  between  completely  and  partially  covered 
situations 

The  comparison  of  these  two  situations  is  shown  in  Fig.  10.  The 
length  of  the  heat  exchanger  (L)  is  fixed  at  180  mm  and  the  length 
of  the  TEG  cuboid  (Lteg)  is  varied  in  the  partially  covered  situation, 
p  Fig.  10  shows  with  the  same  number  of  TE  couples,  the  partial 

g  covered  case  always  performs  better  than  the  completely  covered 

^  one;  the  difference  increases  with  the  increasing  strength  of  the 
heat  sink  and  decreases  with  the  increasing  length  of  the  TEG 
fV  cuboid,  or  with  the  relatively  decreasing  length  of  the  heat 
exchanger.  The  total  power  in  the  partially  covered  case  with 
Lteg  =  80  mm  is  larger  than  that  in  the  completely  covered  case 
with  Lteg  =  180  mm  when  he  =  1800  W/m2  K.  It  implies  a  larger 
power  can  be  generated  with  less  TE  couples  by  simply  making  the 
heat  exchanger  longer  than  the  TEG  cuboid.  Fig.  11  shows  again  that 
the  temperature  difference  between  the  hot  and  cold  sides  of  the 
TEG  cuboid  at  a  fixed  downstream  distance  is  smaller  and  the  up¬ 
trend  near  the  downstream  edge  of  the  TEG  cuboid  becomes  less 
obvious  when  the  TEG  cuboid  is  longer  or  the  uncovered  part  of  the 
heat  exchanger  is  shorter.  All  these  observations  are  consistent 
with  those  found  in  the  last  subsections  and  are  strongly  correlated 
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Fig.  7.  The  distributions  of  the  temperature  difference  between  the  hot  and  cold  sides  of  the  TEG  cuboid  along  the  z  direction  when  the  heat  exchanger  is  partially  covered  with  TE 
couples. 


to  the  fact  that  the  strong  heat  sinks  make  cooler  the  wall  where 
they  are  attached  so  as  to  loot  heat  from  the  neighboring  hotter 
wall.  Consequently,  adding  more  TE  couples  at  the  downstream 
side  generates  an  additional  power  but  also  reduces  the  upstream 
temperature  difference  and  consequently  the  power  generation  of 
the  upstream  TE  couples.  The  additional  power  by  the  downstream 
TE  couples  is  thus  offset  by  the  reduced  power  of  the  upstream  TE 
couples.  The  total  power  gets  saturated  or  even  turns  decreasing 
eventually  as  observed  from  Figs.  6  and  10. 

3.4.  Waste-heat  recovery  efficiency 

Among  all  the  energy  harvesting  systems  performed  herein,  the 
maximum  total  power  (Pmax)  is  about  48  W  and  the  total  heat 
transferred  to  the  TEG  cuboids  is  about  3.3  kW  when  the  heat 
exchanger  of  length  180  mm  is  partially  covered  with  TEGs  of 
length  80  mm  and  the  heat  sink  of  strength  he  =  1800  W/m2  K  is 
employed.  The  input  extractable  energy  may  be  expressed  as 

Qin  =  mcp(Tin-Ta),  (6) 

where  rh  is  the  mass  flow  rate  and  cp  is  the  heat  capacity  of  air  at 
(Tin  +  Ta)/2.  In  all  of  our  simulations,  Q.in  is  8.2  kW.  The  efficiency  of  the 
heat  exchanger  (^exchanger  =  Qabsorb/Qin)  is  about  40%,  the  conversion 


Table  3 

The  total  amount  of  the  heat  transferred  to  the  TEG  cuboid  (Qabsorb)  and  the  part 
transferred  from  the  downstream  side  (Qd). 


L  (mm) 

60 

120 

150 

180 

210 

240 

he  =  600  W/m2 
Qd  (W) 

I< 

95 

226 

243 

250 

253 

250 

Qabsorb  (W) 

1816 

1921 

1933 

1941 

1943 

1938 

Qd/Qabsorb 

5% 

12% 

13% 

13% 

13% 

13% 

he  =  1800  W/m 
Qd  (W) 

2I< 

223 

404 

426 

437 

439 

445 

Qabsorb  (W) 

2354 

2524 

2544 

2559 

2560 

2552 

Qd/Qabsorb 

10% 

16% 

17% 

17% 

17% 

17% 

efficiency  of  the  TEG  (t^teg  =  Pmax/ Qabsorb)  is  about  1.5%,  and  the  re¬ 
covery  efficiency  (^recovery  =  Pmax/Qin)  is  about  0.6%.  These  effi¬ 
ciencies  are  relatively  low  compared  to  the  existing  ones  in  the 
literature,  mainly  because  the  whole  system  is  not  optimized  yet.  For 
instance,  a  replacement  of  the  air  cooling  system  by  a  water  cooling 
system  or  a  use  of  a  highly  efficient  heat  exchanger  (such  as  coun¬ 
tercurrent  heat  exchangers)  can  improve  the  performance  signifi¬ 
cantly.  Nonetheless,  the  present  investigation  highlights  the 
existence  of  an  optimum  coverage  rate  of  TE  couples  on  the  heat 
exchanger. 


Fig.  8.  The  wall  temperature  distributions  of  the  heat  exchanger  along  the  z  direction.  The  section  attached  with  TE  couples  is  colored  by  red.  (For  interpretation  of  the  references  to 
color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  9.  The  total  power  generation  rate  against  the  length  of  the  heat  exchanger  when 
the  length  of  the  TEG  cuboid  (LTEG)  is  fixed  at  50  mm. 


At  last,  we  make  some  discussions  about  the  saving  of  the  en¬ 
gine  shaft  power  by  installing  such  a  TEG  energy  harvesting  system. 
The  installation  of  the  TEG  system  generates  additional  electric 
power,  but  also  causes  some  parasitic  power  loss,  such  as  blow¬ 
down  power  loss  due  to  the  backpressure  increase  of  the  engine 
and  the  power  loss  due  to  the  increased  weight  by  the  TEG  system. 
In  the  maximum-power  case  discussed  above,  the  induced  pressure 
drop  is  about  25  Pa,  much  smaller  than  the  pressure  drop  across  the 
muffler  which  is  about  30  kPa  [24].  Karri  et  al.  [17]  also  indicated 
the  dominant  loss  arises  from  the  increased  weight  of  the  TEG 
system.  According  to  Karri  et  al.  [17],  the  saved  engine  shaft  power 
may  be  evaluated  by  the  following  formula: 


^saved 


^PCU^max  /3rMV 

^alternator  Vd 


(7) 


where  appear  the  rolling  resistance  coefficient  /iT,  the  weight  of  the 
TEG  system  M,  the  speed  of  the  vehicle  V ,  the  driveline  transmission 
efficiency  tjd.  the  Power  Conditioning  Unit  (PCU)  efficiency  ripcu, 
and  the  alternator  efficiency  ^alternator-  By  using  their  parameters  as 
rearranged  in  Table  4  and  substituting  a  vehicle  velocity  of  30  m/s 
and  a  weight  of  4.5  kgw  (the  weight  of  the  maximum-power  sys¬ 
tem  including  the  aluminum  heat  exchanger  and  the  S2090-40W 
heat  sinks  [21])  into  the  formula,  the  saved  engine  shaft  power  is 
about  67  W. 


Fig.  10.  The  total  power  and  the  average  power  per  TE  couple  against  the  length  of  the  TEG  cuboid,  LTEG;  L  =  180  mm  is  fixed  in  the  partially  covered  situation. 


Fig.  11.  The  distributions  of  the  temperature  difference  between  the  hot  and  cold  sides  of  the  TEG  cuboid  along  the  z  direction  when  the  heat  exchanger  length  (L)  is  fixed  at  180  mm. 
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Table  4 

The  values  of  the  parameters  employed  for  evalu¬ 
ating  the  saved  engine  shaft  power. 


Mr  [17] 

0.012 

VD  [17] 

0.9 

Vpoj  [17] 

0.88 

^/alternator  [17] 

0.5 

4.  Conclusion 

In  this  work,  a  heat  exchanger  attached  with  TEGs  for  recovering 
waste  heat  from  an  automotive  exhaust  pipe  was  investigated  by 
using  the  commercial  software  ANSYS-Fluent.  It  was  found  the  heat 
sinks  attached  to  the  downstream  TEGs  possibly  cause  the  down¬ 
stream  wall  cooler  which  thus  loots  heat  from  the  upstream  hotter 
wall,  resulting  in  a  degradation  of  the  performance  of  the  upstream 
TEGs.  Consequently  implementing  more  TE  couples  may  not 
necessarily  be  worthwhile  economically.  Leaving  some  down¬ 
stream  part  of  the  heat  exchanger  uncovered  with  the  TE  couples 
on  the  other  hand  results  in  a  hotter  downstream  wall,  more  heat 
transferred  to  the  upstream  TE  couples,  and  consequently  a  larger 
power  generation  rate.  A  saturation  of  the  power  generation  occurs 
nonetheless  when  the  TEG-uncovered  part  of  the  heat  exchanger  is 
sufficiently  long.  This  investigation  reveals  a  new  parameter  to 
optimize  the  energy  harvesting  system— the  number  of  the 
employed  TEG  chips,  once  the  heat  exchanger,  the  cooling  system, 
and  the  thermoelectric  material  have  been  carefully  determined. 
For  the  heat  exchanger  investigated  herein,  the  optimum  length  of 
the  TEG  cuboid  is  80  mm  for  a  heat  exchanger  of  length  180  mm, 
and  the  optimum  length  of  the  heat  exchanger  is  120  mm  for  a  TEG 
cuboid  of  length  50  mm,  when  the  convection  heat  transfer  coef¬ 
ficient  associated  with  the  heat  sink  is  1800  W/m2  K. 
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